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Abstract 

We follow the approach used by Davies and Pringie (1981) and discuss the 
trans-sonic substage of the propeller regime. This substage is intermediate 
between the supersonic and subsonic propeller substages. In the trans-sonic 
regime an envelope around a magnetosphere of a neutron star passes through 
a kind of a reorganization process. The envelope in this regime consists of 
two parts. In the bottom one turbulent motions are subsonic. Then at some 
distance rs the turbulent velocity becomes equal to the sound velocity. During 
this substage the boundary rs propagates outwards till it reaches the outer 
boundary, and so the subsonic regime starts. 

We found that the trans-sonic substage is unstable, so the transition be- 
tween supersonic and subsonic substages proceeds on the dynamical time scale. 
For realistic parameters this time is in the range from weeks to years. 
Keywords: compact objects, isolated neutron stars, evolution 

1 Introduction 

Observational appearances of a neutron star (NS) are mainly determined by its 
interaction with the surrounding plasma. The following main stages (regimes) can 
be distingushed (see a very detailed description in Lipunov 1992 or in Lipunov et 
al. 1996) 



• Ejector. At this stage plasma is swept away by low-frequency electromagnetic 
radiation or/and by a flow of rclativistic particles. Matter is stopped further 
than the so-called light cylinder radius rg. 

• Propeller. If a NS is in the propeller regime than matter can penetrate inside 
re, but it is stopped by a rapidly rotating magnetosphere of the NS. 

• Accretion. Finally the NS is slowed down and the centrifugal barrier disap- 
peares, so if matter cools fast enough then it can fall down onto the surface 
of the compact object. 



Normally a NS is born at the stage of ejection (a radio pulsar is a classical ex- 
ample of a NS at the ejector stage). Then as the spin period increases the NS passes 
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Figure 1: Pressure vs. radius. Left panel: external pressure. Right panel: internal 
pressure. 

propeller and accretor stages. For NSs with large (>400 km/s) spatial velocities 
another stage — Georotator — can appear. 

In a simplified model it is possible to define transitions between different stages 
by comparing external and internal pressure (Fig. The external one can be 
roughly approximated as the ram pressure of a flow of the interstellar medium 
(far from the NS) or as the pressure of matter falling down onto the NS in its 
gravitational field (for distances smaller than the so-called gravitational capture 
radius tq). The internal one inside the light cylinder ri can be estimated as a 
pressure connected with the magnetic field of a NS. 

In this paper we discuss the propeller stage. An existence of it was recog- 
nized long ago (initially by Shvartsman (1970) and later on by lUarionov and Sun- 
yaev(1975). However still this stage is not well understood. Here we introduce and 
discuss a new substage of this regime of magneto-rotational evolution of NSs. 

2 General features of the propeller stage 

Davies et al. (1979) and Davies, Pringle (1981) note that at the stage of propeller 
there can be a significant energy release at the magnetospheric boundary. Energy 
can be large enough to form a kind of a turbulent quasistatic atmosphere. These 
authors distinguished three substages of the propeller regime. 

1. Very rapid rotator: 

Cs(nn) ^ Ti^Q. > Wff . (1) 

here Cg - sound velocity, = 27r/P - spin frequency, ve - free-fall velocity. 



2. Supersonic propeller: 



ri„r2 > Cs(nn) • 



(2) 
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3. Subsonic propeller: 

rinO < Cs(rin) (3) 

and 

Vt{r)<Cs{r) in nn < r < rout = m , (4) 

here and rout — are internal and external radii of an envelope, Vt — turbulent 
velocity. 

We will not discuss the substage of the very rapid propeller here. Mainly we 
focus on super- and subsonic substages and on the transition between them. The 
supersonic regime can be considered a classical propeller where accretion is impossi- 
ble due to a centrifugal barrier. At the subsonic stage the magnetospheric (Alfven) 
radius is smaller than the corotation radius re- Accretion does not start because 
temperature is too high (see original discussion in Davies, Pringle, 1981 and later 
proposals in Ikhsanov 2003). 

Braking laws (for the spin-down) are different at different substages. In the su- 
personic regime energy loss rate does not depend on the spin period (Davies, Pringle 
1981), however different formulas for the spin-down at this stage were suggested, 
see a list for example in (Lipunov, Popov 1995). In the subsonic one the spin-down 
is always slower: P <xt. 

3 Why does the intermediate regime exist? 

In general supersonic and subsonic regimes cover all possible values of the rotation 
velocity = 27r/P. The supersonic propeller formally operates till rl^^'^^Vl > 

Cs(C"') • Here rf^"" = r^J^'rH^ (Davies, Pringle 1981), ru = {ii'' / {MV2GM)f/^ 
- magnetospheric (Alfven) radius. 
The subsonic regime is on when 

C'f^<c,(rf„"'') and 7.t(r) <Cs(r) for rf„"^ < r < rout - re , (5) 

here rf^'' = ryi (Davies, Pringle 1981). As for the subsonic stage vt{rf^'') ~ rfj^'Ti 
and Wt/Cs oc r^/^ , then this regime is valid for rfj^'Ti < Cs(rf^"'') {rf^^ /ro)^^^ . 
Note the following properties of the stages. 

1. At both stages internal radii of an envelope rin do not depend on O, and 
always rf„"^'='' > rf^K 

2. The structure of the envelope on the two substages is different. 

3. It is easy to check that the end of the supersonic substage and the beginning 
of the subsonic one both correspond to: 

O = V2GMr^'\^^'\ (6) 
P-l.lSlO^s/i^/VVlaf- 
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Even in the frame of Davies and Pringle model an intermediate regime during 
which the structure of the envelope is reorganized is inevitable. We call this interme- 
diate stage trans-sonic propeller. As wc show below this is a short noncquilibrium 
episode in a life of a NS, and the spin frequency does not change significantly during 
this reorganization process. 



4 Trans-sonic propeller 

Let us consider the structure of a quasistationary envelope at the intermediate trans- 
sonic propeller substage. In the atmosphere around such a NS Cs(r) ~ vg (Davies, 
Pringle 1981). (This condition is valid also for super and subsonic propeller.) 
Processes in the lower part of the atmosphere are similar to the ones on the subsonic 
stage: 

Vt{rin) ^ Tinfl < Cg. (7) 

If the envelope (atmosphere) is adiabatic then for this bottom part of it the 
politropic index is equal to n = 3/2 and 

p(r) oc ^ p(^) ^ J.-5/2 ^ (-g^ 

We assume following Davies and Pringle (1981) that the rotational energy of the NS 
is dissipated at the magnetospheric boundary and that it is transported outwards 
by turbulence. For such assumptions we have: 

vt{r) oc r-i/*' (9) 

and the turbulent Mach number is: 

M(r) = ^ocrV3. (10) 

Till Ait < 1; i-e. while r < Tg where Tg is the boundary between two parts of 
the envelope (rin < < tg) the structure of the envelope is not changed. For large 
radii turbulence becomes supersonic. Small-scale shock waves are formed and they 
quickly dissipate part of the energy, so that turbulent velocity decreases down to 
the sound velocity. In the range rg < r < re the envelope structure is different from 
the bottom part: 

Xt(r)~l, (11) 
p{r) oc , p{r) oc r'^/^ . (12) 

In the outer part of the envelope physical conditions are similar to the ones on 

the supersonic substage. 

To determine parameters of the whole atmosphere it is necessary to calculate 
the position of the boundary between two parts of the envelope, Vg, and position of 
the inner boundary of the bottom part, (during the transition it decreases from 
^G^^M^ to tm). To do it it is necessary to solve the following system of equations: 
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o 7/6 -1/6 I2GM 



— 7/2 

However, the system is degenerate and each equation can be reduced to: rg (x rj^^ . 
If the following equation is fulfilled: 

^2 / 1 Mv^ \ (2GAf)3/2 



(14) 



then the system is compatible, i.e. for all in the range ryi < rjn < ^m^i^q^ 
There is some Vs that is a solution of the eq. ^| 

The compatibility condition is fulfilled at the end of the supersonic substage. 
Later (during the transition) at any given moment (for any Q) left-hand side of eg 1141 
is smaller than the right-hand one. It means that the magnetospheric pressure and 
the envelope pressure have the same dependences on rjn, but the latter one is always 
larger (the first equation of the system ^J. 

During the trans-sonic stage the period is not changing significantly (a typical 
value is determined by eq. 6), so in terms of the rotational evolution the subsonic 
substage nearly immediatclly follows the supersonic one. The spin-down law for the 
trans-sonic propeller is the same as for the subsonic regime. 

An energy release during the transition stage is negligible. Estimates for realistic 
isolated NS parameters give a value AE < 10^° erg. 



5 Discussion 

We want to note that calculation similar to the ones presented above are just rough 
estimates. There are several reasons for that. 

The first is connected with uncertainties in many parameters, even in their 
determinations. For example the accretion rate is usually taken as M = TrrQPooVoc- 
This is just an estimate, and for different velocitites it can be different from the 
actual value by a factor a few. Small changes is some parameters can lead to 
significant changes in others. For example Ikhsanov (2001) discusses the value of 
the critical period which determines the end of the subsonic propeller stage (and so 
the accretion stage begins). The obtained value is different from the one found in 
the original paper by Davies and Pringle (1981) by a factor of 7.5. Correspondently 
all time scales are also significantly changed. But note, that this fact is due to a 
change in the magnetospheric radius {tm) only by a factor of '--^ 2! 

The second one is connected with idealizations. Even if all parameters can be 
well defined, then it is necessary to take into account such details as non-spherical 
form of the magnitosphere, inclination of the magnetic axis relative to the spin 
axis, angular moment in the infalling matter (even for cases when the condition 
for disc formation is not fulfilled), etc. For example even if t^m > t'c part of the 
magnetosphere is inside Tc as the corotation radius is the radius of a cylinder, not a 
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sphere. In that sense the process of ahgnment (see for example Regimbau, de Freitas 
Pacheco 2001 and Beskin et ah 1993) during the magneto-rotational evolution of a 
NS can be important in a destiny of a NS. 

A low rate accretion can proceed even at the propeller stage due to several 
reasons. One of them is diffusion of plasma. Such form of accretion was discussed in 
details by Ikhsanov (2003). For long spin periods luminosity due to such an accretion 
can be larger than the dissipation of the the rotational energy on the boundary of 
the magnetosphere. 

An important question is connected with the whole time of evolution prior the 
accretor stage ^a- Obviously 

tA = tE + tp, (15) 

where is the time which a NS spends at the stage of ejector, and tp is the 
duration of the propeller stage. Even is is not well determined. Usually authors 
assume that spin-down at this stage is determined by the magneto-dipole formula 
with the braking index equal to 3. However, direct measurement for many radio 
pulsars show that the braking index is smaller than 3. Also an evolution of the 
angle between spin and magnetic axes is usually not taken into account. 

As the propeller stage consists of the three substages then it is necessary to 
write: 

^supcr ^" ^trans ^ ^sub- (-i^) 

If t-E and tp are determined the fate of NSs for different parameters can be 
easily shown on — tp-diagram suggested in (Popov 2004). In the FigI21we show 
an example of such a plot. For this illustration we assume that the accretion regime 
starts at the period: 

Pb. = 4.5 . 10Vj^/^^M-^/^m-^/^\ (17) 

Ms = Af/10^gs~^, m — My\'s,IMq. Here we renormalize the value by Ikhsanov 
(2001). 

Time scales are determined by (see Popov 2004): 

tE = 0.8-10V3o'""'^'«6yr, (18) 
here n = pm~^ - interstellar medium number density, nip - proton mass. 

t,upcr = 1.3 • lO^fi^o'^'n-^/'vl^' yr, (19) 

this is a very efficient spin-down suggested by Shakura (1975). so our estimate of 
4upcr is in some sense a low bound. 

itrans + ^sub = 10^ ^J■3o'm Phr yr. (20) 

The mass of a NS is assumed to be Mns = 1.4 Mq. 

In the plot we show lines for ^e + tp equal to 1, 5 and 10 Gyrs. Eight symbols 
corresponds to eight combitations of n,v,fj. (see the table). 



Trans-sonic propeller stage 



7 



11 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 


1 1 1 1 1 1 1 1 1 1 


10 




+ 7 - 






o 


>> 

J- 9 


't 


+ \ 

7 \ 




Log 


+4 






8 


- 

+ 6 i 






7 


1 ' 


i 1 1 






7 8 9 10 11 

Log tg, yrs 



Figure 2: — fp-diagram for isolated NSs. Quasi-linear nature of the distribution 
of model points on the log ts - log tp plane was discussed in detail in (Popov 2004). 

6 Conclusions 

In this paper we showed that an intermediate substage of the propeller regime - 
the trans-sonic propeller - should exist. However, the stage is non-stationary and 
very short. As it was shown above the existenee of this intermediate stage does not 
change the timescale of the evolution prior to the accretor stage (see Ikhsanov 2003 
for a discussion of the timescales). Other conclusions can be summarized as follows: 

• The intermediate trans-sonic propeller substage in unstable. 

• The duration of the transition can be roughly estimated as ^rc/vff (from 
weeks to years for realistic isolated NSs). 

• The spin frequency is nearly unchanged during this transition. 

• The energy release during the transition is small. 
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Table 1: t-E and tp for typical values ofn, v and fi 



Number n v, km s 



^30 log tE, yrs log tp, yrs 



1 
2 
3 
4 
5 
6 
7 
8 



0.1 
0.1 
0.1 
0.1 
1.0 
1.0 
1.0 
1.0 



20 
20 
40 
40 
20 
20 
40 
40 



1 

10 
1 

10 
1 

10 
1 

10 



9.70 
8.70 
10.0 
9.01 
9.20 
8.20 
9.51 
8.51 



9.39 
8.15 
10.0 
8.79 
8.68 
7.44 
9.32 
8.08 
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